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Abstract 

Based on 58 x 10 6 J/ip events collected by BESII, the decay J/ip — > K^Ksit^it is 
studied. In the invariant mass spectrum recoiling against the charged K* (892) 1 * 1 , the 
charged n particle is found as a low mass enhancement. If a Breit-Wigner function 
of constant width is used to parameterize the k, its pole locates at (849 ± 771^) — 
z(256 ±40i$j) MeV/c 2 . Also in this channel, the decay J/tp -> iT* (892) + IT* (892) ~ 
is observed for the first time. Its branching ratio is (1.00 ± 0.19^g 32) x 10~ 3 . 



1 Introduction 

The a and k are controversial particles in hadron spectroscopy. They were first 
found in the analysis of nir and ttK scattering data. Because the total phase 
shifts in the lower mass region are much less than 180° and they do not fit 
into ordinary qq meson nonets, they have been the subject of violent debates. 
Refs. ([T] - [12]) are some recent analyses that support their existence. 

Evidence for k particles comes from the study of production processes and the 
re-analysis of Kit scattering data. Evidence for the k has been found in D + — >■ 
K-TT+7T- |BJ, J/ij) i?"*(892)°ii" + 7r- 0E], and D+ K~^^v [13]. A Kir 
s-wave component is found in D° — > K~K + 7r° [H], D + — > K~ir + e + i> e |15j . 
and t — > Ks7r~u T [16]. But no evidence for the k in D° — > K~7t + tt° [17] 
or for the charged k in D° — > K~K + tt° [18j is seen. The k is found in the 
phenomenological analysis of Kit scattering phase shift data (p2], [19] - |26|). 
However, some theorists are not convinced by the evidence ([27] - [30]). The 
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present status of the k is summarized by the Particle Data Group PDG [3T] . 

The a and k particles have been studied with BES data, where evidence for 
a and k particles is quite clear ([1] - [6]). A neutral k was found in the decay 
J/ip -» K*(892)°k° -» K*(892) K + %- 0(6]. Because of isospin symmetry, if 
a neutral k exists, a charged n should exist and could be produced in J/ip — > 
K*(892) ± k^. 

In this study, we search for and study the charged k in J/ip — > K Ks^7^° ■ 
Our analysis is based on 58 million J/ip decays collected by BESII at the 
BEPC (Beijing Electron Positron Collider). BESII is a large solid-angle mag- 
netic spectrometer which is described in detail in Ref. [32]. The momentum 
of charged particles is determined by a 40-layer cylindrical main drift cham- 
ber (MDC). Particle identification is accomplished using specific ionization 
(dE/dx) measurements in the MDC and time-of-flight (TOF) information in 
a barrel-like array of 48 scintillation counters. Outside of the TOF is a barrel 
shower counter (BSC) which measures the energy and direction of photons. 



2 Event selection 

In the event selection, candidate tracks are required to have a good track fit 
with the point of closest approach of the track to the beam axis being within 
the interaction region of 2 cm in r xy and ±20 cm in z (the beam direction), 
polar angles 9 satisfying |cos#| < 0.80, and transverse momenta Pt > 50 
MeV/c. Photons are required to be isolated from charged tracks, to come 
from the interaction region, and have deposited energy in the BSC greater 
than 40 MeV. Events are required to have four good charged tracks with total 
charge zero and at least two good photons. 

For the K$ reconstruction, we loop over all oppositely charged pairs of tracks, 
assuming them to be pions, and fit them to K$ —> 7r + 7r~, which determines 
vertices and 7r + 7r~ invariant masses, M n7T , for the four possible combinations. 
The combination with M n7T closest to Mk s is selected and is required to satisfy 
\M n + n - — M Ks \ < 20 MeV/c 2 and have its decay vertex in the xy-plane satisfy 
r xy > 0.008 m. After K$ selection, the particle type of the remaining two 
tracks, that is whether they are K + tt~ or K~tt + , is decided by selecting the 
combination with smallest Xtof + Xdedx- 

A four constraint (4C) kinematic fit is applied under the i ; T ± 7r T 7r + 7r~77 hy- 
pothesis, and xlc < 20 is required. Events with a 77 invariant mass satisfying 
|M 77 -M n o \ < 40 MeV/c 2 are fitted with a 5C kinematic fit to ^Tr^n+n-n 
with the two photons constrained to the ir° mass, and events with \\c < 50 
are selected. The 7r T 7r° mass distribution is shown in Fig. 1(a), where the 
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p is clearly seen. Decays with an intermediate p are background, and the 
requirement \M n ^ n o — M p \ > 100 MeV/c 2 is applied to remove them. The 
requirements \M Ks7T o - 0.897| > 40 MeV/c 2 and \M K ± wT - 0.897| > 40 
MeV/c 2 are used to remove backgrounds from J/ip — > K*(892)°K ± ir T and 
J/tp — > K* (892)° K sir . After these requirements, the combined K ± n° and 
Ksir T mass distribution is shown in Fig. 1(b); the highest narrow peak is 
charged If* (892). The M K ± n o versus M Ks7T ^ scatter plot is shown in Fig. 1(c). 
There are two clear bands, a vertical and horizontal band, which correspond 
to J/ip K*(892) ± K s n T and J/ip -> if* (892) ±7^77°, respectively. The re- 
quirements \M K±7T o - 0.892| < 80 MeV/c 2 and \M Ks7T ± - 0.892| < 80 MeV/c 2 
are imposed to select J /if; ->• K*(892) ± K s n^ and J/ip K*(W2) ± K*n a 
events, respectively. 

After the above selection, the combined Ksn T and K t tt° invariant mass dis- 
tribution recoiling against the K* (892) ± is shown in Fig. 1(d). It is the sum of 
the four decays listed in Table 1 with about 1000 events for each, as listed in 
the table, and a total of 4121 events. In Fig. 1(d), a clear narrow peak at 892 
MeV/c 2 and a wider peak at about 1430 MeV/c 2 are seen. In addition, there 
is a broad low mass enhancement just above threshold. The spectrum is quite 
similar to the spectrum of K + tt~ in the decay J/ip — > K*(892)°K + tt~ [6]. The 
biggest difference between them is that the charged i^*(892) peak is much 
larger than the neutral i^*(892) peak. The i^*(892)7r invariant mass distribu- 
tion is shown in Fig. 1(e), where there are indications of peaks at 1270 MeV/c 2 
and 1400 MeV/c 2 . The resulting Dalitz plot is shown in Fig. 1(f). The two 
diagonal bands correspond to the low mass enhancement combined with the 
892 MeV/c 2 peak and the peak around 1430 MeV/c 2 in the Kit spectrum. 



Channel 


Events 


J/ip -> K* (892)+ K s -K~ -> K + tt Kstt- 


1023 


J/ip -> K*(892)-K s ir + -> K--K K S Tr + 


946 


J/ip -> K*(892)+K-tt° -> Ks-ir+K-iT 


1055 


J/ip -> K*(892)-i^+vr° -)■ K s tt-K + -k° 


1097 



Table 1 

Four signal channels and number of events from each. 



3 Background Studies 

Possible sources of background are studied. First, sideband backgrounds are 
studied. The 7r + 7r~ mass distribution is shown in Fig. 2(a), where a clear K $ 
signal can be seen and the background level is quite low. The Kir spectrum 
from Ks side-band events is shown in Fig. 2(b). (The K$ side-band is defined 
by 0.02 GeV/c 2 < |M W7r -0.497| < 0.04 GeV/c 2 .) There are no clear structures. 
The 77 spectrum is shown in Fig. 2(c), and the Kn spectrum of the 284 7r° side- 



4 



N jooo 



> 
a> 



^800 
o 

CM 

1 600 

> 
hi 



400 



200 




8.20.40.60.8 1 1.21.41.61.8 2 

(GeV/c 2 ) 




.60.8 1 1.21.41.61.8 2 2.2 

M K n (GeV /c 2 ) 



> 
a> 
o 



2.4 
2.2 

2 
1.8 
1.6 
1.4 
1.2 

1 

0.8 
0.6 





(c) 








i 




:-S'n;::">. : :^:i;-:v :t 















0.60.8 1 1.21.41.61.8 2 2.22.4 

(GeV / c 2 ) 



M 




K+- Jt0 



0.60.8 1 1.21.41.61.8 2 2.2 

M Kjt (GeV/c 2 ) 




1 1.21.41.61.8 2 2.22.42.6 

M,p„ (GeV/c 2 ) 



o 

> 
d> 
O 



CO 




2 3 4 5 6 7 8 9 

M K*(892) K ( QeV/C ) 



Fig. 1. (a) tt^tt mass distribution, (b) Combined K^ir and Xstt 1 ' 1 mass distri- 
bution after final cuts except the iT*(892) cuts, (c) Scatter plot of M K ± n o ver- 
sus M Ks7T ^. (d) Invariant mass distribution of K^ir and recoiling against 
K*(892) =F . (e) Invariant mass distribution of K*(892)vr. (f) Dalitz plot. 
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band events is shown in Fig. 2(d). (The n sideband is defined by 0.04 GeV/c 2 
< |m 77 — 0.135| < 0.08 GeV/c 2 .) The structures in Fig. 2(d) are similar to those 
in the signal region and come from signal events in the ir° tails. The K* (892) 
side-band background is shown by the dark shaded histogram in Fig. 2(e). 
(The K*(892) side-band is defined by 0.08 GeV/c 2 < \M K7T - 0.892| < 0.16 
GeV/c 2 .) The clear K*(892) in the Kn mass distribution of K*(892) side- 
band events mainly comes from the cross channel. (There are two bands in 
the scatter plot in Fig. 1(c). When we select one band, we will also select 
some events from the other band where it crosses the first band. These events 
correspond to cross channel background.) Fig. 2(f) shows the Kir spectrum 
after side-band subtraction. The low mass enhancement and K*(892) peak 
survive after side-band subtraction. 

Next, we perform Monte Carlo simulation to study the main physics back- 
ground processes, including J/ip — > jr/c — > jK*K* — > •yK ± Ks'n T iT , J/ip — > 
77/ c 7^(892)^^° ->- jK+Ksn^ir , J/ip jr] c ->• 7^(892)^7^ 
— > / ~fK ± KsTi' T 7T , J/ip — > 7r°7r + 7r~7r + 7r - , and J/ip — > 7r K + K~7T + 7T~ . The se- 
lection efficiencies are much lower than 1%, and the largest number of back- 
ground events contributed is about 6. Therefore, the physics background is 
quite low. 



4 Partial wave analysis 



A partial wave analysis (PWA), which is based on the covariant helicity am- 
plitude analysis ([33] - [31]), i s performed for the charged n. We add the 
likelihoods of all four channels together, and find the minimum of the sum. 
This is the same method as used to fit J/ip — > i^*(892) i^ + 7r _ [6]. The main 
difference here is that the decay J/ip — > K*(892) ± K*(892) T is included in the 
fit. 

In the PWA analysis, ten resonances, k, K*(U30), IPS, K*{U30), K*{1920), 
if* (1410) and K*(892) in the Ktt spectrum, 7^(1270) and 7^(1400) in the 
ir*(892)7r spectrum, and &i(1235) in the i^*(892)i^ spectrum, and two back- 
grounds, listed as the last items in Table 2, are considered in the fit. In Table 2, 
IPS (Interference Phase Space) refers to a broad + structure with a Ktc in- 
variant mass spectrum the same as phase space, that interferes with k. K* BG 
refers to the K*(892) background coming from the cross channel. PS (Phase 
Space) refers to the background with no interference with resonances and with 
a shape almost the same as that of phase space. 
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Fig. 2. (a) m w + w - mass distribution after final event selection except for the K$ 
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selection except the 7r° requirement, (d) Mk-k distribution from ir° side-band events, 
(e) Kir spectrum recoiling against if* (892). The dark shaded histogram is from K* 
side-band events, (f) The Kir spectrum after side-band subtraction. 
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Resonance 


Spin- 


Decay 


Mass 


Width 


Sig. 




Parity 


Mode 


(MeV//cr) 


(MeV/cf ) 




K (1) 


+ 


Kir 


810 ± 68^24 


536 ± 87l47 6 


> 6(7 


K (2) 


+ 


Kir 


884 ± 40^22 


478 d= 77„4^ 


> 6(7 


K (3) 


+ 


Kir 


1165 ± 58l4i° 


1349 ± 500+^g 


> 6(7 


i^o (1430) 


+ 


Kir 


1400 ± 86 


325 ± 200 


0.6(7 


IPS 


+ 


Kir 






> 6(7 


#1(1430) 


2 + 


Kir 


1411 ± 30 


111 ± 46 


> 6(7 


if* (1920) 


2+ 


Kir 


2020 ± 140 


705 ±160 


> 6(7 


if*(1410) 


1~ 


Kir 


1420 ± 14 


130 ± 28 


> 6(7 


if* (892) 


1~ 


Kir 


896 db 8 


57 db 12 


> 6(7 


Ki (1270) 


1^ 


K (892)7r 


"1 O f A 1 1 A 

1254 ± 14 


60 ± 28 


> 6(7 


if 1 (1400) 


1 + 


if*(892)vr 


1390 ± 30 


146 ± 44 


> 6(7 


6i(1235) 


1 + 


if*(892)if 


1230 ± 52 


142 ± 38 


4.5(7 


if* (892) BG 


l - 


Kir 






> 6(7 


PS BG 













Table 2 

Resonances included in the fit of this channel. Masses and widths of various reso- 
nances are determined by mass and width scans, k (1), (2) and (3) are results given 
by fits using Breit-Wigner functions (1), (2) and (3) to fit k respectively. IPS refers 
to the broad + structure which interferes with k. K* BG refers to the if* (892) 
background coming from the cross channel. PS BG refers to the background with 
no interference with resonances. 

Three different parameterizations are used to fit the k. They are 
1 



BW K 
BW R 
BW K 



. -— , T K = constant, 
mi — s — imA K 



ml- s - iy/sT K (s] 



r*(s) 



9l-K 

Sirs 



m\ - s - i^T K (s) 



TJs) = a-k K , 



(1) 
(2) 
(3) 



where k K is the magnitude of the K momentum in the Kir, or the k, center 
of mass system [38J, and a is a constant which will be determined by the fit. 
Parameters in the Breit-Wigner function are determined by mass and width 
scans. The minima of the scan curves give the central values of mass and width 
parameters. From these, the corresponding Breit-Wigner pole positions can be 
directly calculated from equation (1), (2) and (3). Our final results are listed 
in Table 3, where the first errors are statistical, and the second are systematic. 
The mass and width parameters obtained by different parameterizations are 
quite different, but their poles are almost the same, which is quite similar to 
what was found in the study of the neutral k. The results for the neutral k [6] 
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are shown in Table 4 and are consistent with the charged k. 





BW (1) 


BW (2) 


BW (3) 


Mass (MeV/c 2 ) 


810 ± 68i^ 


884 ± 40tH 


1165 ± 58±lf° 


Width (MeV/c 2 ) 


536 ± 87+_l^ 


478 ± 77+l\ 


1349 ± 500l^g 


pole (MeV/c 2 ) 


(849±77l^) 
-i(256 ± 40±f|) 


(849 ± 511^) 
-t(288 ± 1011^) 


(839 ± 1451 24 ) 
— i(297±51±fg) 



Table 3 

Masses, widths and pole positions of the charged k. In the table, the first errors are 
statistical, and the second are systematic. BW (1) means equation (1) is used to fit 
the k. BW (2) and BW (3) have similar meanings. 





BW (1) 


BW (2) 


BW (3) 


Mass (MeV/c 2 ) 


745 ± 26+1^ 


874 ± 25j^ 


1140 ± 39j^ 


Width (MeV/c 2 ) 


622 ± 77l^g 


518 ± 65±|? 


1370 ± 156±f21 


pole (MeV/c 2 ) 


(799 ±37t$) 
-i(290 ± 33+H) 


(836 ± 38t$) 
-i(329 ± 66i^) 


(811±74± 83 7 ) 
-i(285 ± 20+H) 



Table 4 

Masses, widths and pole positions of the neutral k [6]. BW (1) means equation (1) 
is used to fit the k. BW (2) and BW (3) have similar meaning. 

Our final results correspond to the solution with the minimum least likelihood. 
Differences among solutions with similar likelihood values are included as sys- 
tematic uncertainties. Also included are the effect of removing iT (1430), IPS, 
&i(1235) and IT* (892) IT* (892), the result of a fit with the If* (892) background 
level floating, and the result from a fit using direct side-band subtraction. 

The masses and widths of all resonances obtained by mass and width scans 
are shown in Table 2. In the fit, the contribution from 7^(1430) is small; its 
statistical significance is only 0.6a. Because it is expected in this channel, it is 
included in the final solution. The Kir mass distribution is shown in Fig. 3(a), 
where points with error bars are data, and the light shaded histogram is the 
final fit. In the figure, the dark shaded histogram shows the contribution of 
the charged k. Fig. 3(b) shows the fit for the i^*(892)7r spectrum, and Fig. 4 
shows the angular distributions. 



5 Branching ratio measurements 



The decay J/ip — > K*(892)+k+ contributes 655 events. Monte Carlo simula- 
tion of J /if) -> K*(892)+k+ K+Ksn+ir determines an efficiency of 2.33%, 
and the branching ratio of J/ip -» K*(892)+k~ or J/ip -)■ K*(892)~k+ is 

M= 2.33% xsTxVx 1/9 = ' 1 - 09±ft " > ^» Xl °" (4) 
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Fig. 3. (a) Final fit results for the Kn spectrum. Points with error bars are data, 
the light shaded histogram is the final global fit, and the dark shaded histogram 
is the contribution of the k. (b) Final fit of iC*(892)7r spectrum. Dots with error 
bars are data, and the histogram is the final global fit. There are two peaks in the 
lower mass region. The lower one is fit by the iTi(1270), and the higher one is by 
the If i (1400). 

where, the first error is statistical, the second error is systematic, the factor 
| is because the events are the sum of four channels, and | is the isospin 
factor. In the previous study of the decay Jlip — > K* (892)° K + it~ , the corre- 
sponding branching ratio for the neutral ft Li] is (0.52 - 0.97) x lO" 3 . The two 
results are consistent with isospin symmetry. The systematic error includes un- 
certainties from multi-solutions, from different background fit methods, and 
from removing some components from the fit (^(1430), IPS, &i(1235), and 
J fsT*(892) ± J fs:*(892) T ). 

Also interesting is the existence of the J /if) electromagnetic decay J /if) — > 
K*(892) ± K*(892) T . The peak at 892 MeV/c 2 in the Kn invariant mass spec- 
trum (see Fig. 1(d)) comes from both cross channel background and from 
J/if> — > K*(892) ± K*(892) T . The number of events from the cross channel can 
be calculated approximately from the if* (892) side-band structure, shown in 
Fig. 2(e), where the narrow peak at 892 MeV/c 2 is clear. In the final fit, the 
decay J /if) i^* (892) ± iT* (892) ^ contributes 323 events. The Monte Carlo 
simulation of the decay J/tf) if *(892) ± £f*(892) T ^Ks^n yields an 
efficiency of 1.25%, and the branching ratio of J/ip -> K* (892)+ K* (892)~ is 



BR 



323/2 



1.25% x 5.8 x 10 7 x 1/9 x 2 



'1.00 ±0.19 



+0.1U 

0.32J 



x 10 



-3 



(5) 



1 The branching ratio of the neutral k was not reported in Ref . [6] . In the study of 
the neutral k, the number of k events was in the range 1891 - 3516, and the selection 
efficiency was 14.2%. Its branching ratio is 14 2%x5 1 8xiq 7 x4/9 



(0.52-0.97) x 10- 
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Fig. 4. Final fit result for the angular distributions. Dots with error bars are data, 
and the histogram is the final global fit. The upper left figure shows the final fit of 
the 6\ distribution, the upper right shows the final fit of the <p\ distribution, the 
lower left shows the final fit of the 62 distribution, the lower right shows the final 
fit of the 4>2 distribution. Q\ and <p\ are the polar angle and azimuthal angle of the 
Kit system in the J ftp center of mass system. 62 and (f>2 are the polar angle and 
azimuthal angle of the K meson in the Ktt center of mass system. 



where, the first error is statistical, the second error is systematic, the factor 
I is because events are counted twice, | is the isospin factor, and 2 is be- 
cause the data comes from two decay channels: J/ip — >■ i^*(892) + f^*(892)~ — > 
(K + iT )(K s ir-) and J/ip ->• K* (892)- K* (892)+ (K-tt°)(K s tt + ). The sys- 
tematic error includes uncertainties from multi-solutions, from different back- 
ground fit methods, and from removing some components from fit (^(1430), 
IPS, and 6i(1235)). 
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6 Summary 



In conclusion, the charged k is observed and studied in the decay J '/ip — >■ 
K*(892) ± k t — > K^KsttTtt ■ The low mass enhancement in the Kir spectrum 
cannot be fit well unless a charged k is added into the solution. If we use a 
Breit-Wigner function of constant width to parameterize the k, its pole locates 
at (849 ± 771J!) -i(256 ±4012!) MeV/c 2 . In our analysis, three different « pa- 
rameterizations are tried in the fit, and final results are shown in Table 3 and 
are consistent with those of the neutral k and are also in good agreement with 
those obtained in the analysis of Kir scattering phase shifts. Also, the decay 
J/ip — > i^*(892) ± i ; r*(892) =F is observed for the first time with the branching 
ratio (1.00±0.19±^) x 10 3 . The corresponding decay mode is not observed 
in J/ip -> K*(892)°K + 7i-. The decays J /if) -> iT*(892) ± K*(892)T can be pro- 
duced through J/ip electromagnetic decays, while J/ip -> K* (892)° K* (892)° 
can only be produced through J/ip hadronic decays, which would be SU(3) 
symmetry breaking decays and are suppressed. 
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